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Abstract 


This paper studies the decay of a large, closed domain wall in a closed 
universe. Such walls can form in the presence of a broken, discrete symmetry. 
We study a novel process of quantum decay for such a wall, in which the 
vacuum fluctuates from one discrete state to another throughout one half of 
the universe, so that the wall decays into pure field energy. Equivalently, the 
fluctuation can be thought of as the nucleation of a second closed domain wall 
of zero size, followed by its growth by quantum tunnelling and its collision 
with the first wall, annihilating both. We therefore study the 2-wall system 
coupled to a spherically symmetric gravitational field. We derive a simple form 
of the 2-wall action, use Dirac quantization, obtain the 2-wall wave function 
for annihilation, find from it the barrier factor for this quantum tunneling, 
and thereby get the decay probability. This is the second paper of a series. 
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I. INTRODUCTION 


It is well appreciated that domain walls can wreck a cosmological model, and therefore it 
is of interest to hnd processes that can destroy domain walls in the early universe. In paper 
I [H we introduced a novel process, the quantum decay of domain walls by global fluctuation 
and quantum tunnelling. We studied a closed universe dominated by a single closed domain 
wall — the Vilenkin-Ipser-Sikivie or “VIS” solution 1^,^ — and we found an instanton 
that mediates its decay into a closed universe containing pure held energy. However some 
technical problems cropped up in the instanton calculation. Therefore, in this paper we will 
study the same decay process by a different technique, namely a Hamiltonian formulation 
and Dirac quantization. 

We must hrst explain why gravity is involved in this decay at all. A domain wall in hat 
spacetime separates two inhnite regions of diherent discrete vacuum state, in the presence 
of a broken discrete symmetry. The wall cannot decay because any quantum huctuation 
into a no-wall state has an inhnite barrier. The Vilenkin solution seems to describe 
an inhnite domain wall dressed by its gravitational held; however this spacetime is not 
geodesically complete, and its complete analytic continuation [Q can be interpreted as a 
closed, topologically universe dominated by a closed hnite domain wall. This universe 
(we call it the VIS solution) starts at inhnite volume, collapses to a minimum volume, at 
which point it halts and then re-expands to inhnite volume. The minimum radius of the 
domain wall is Rmin ~ l/crG where a is wall surface tension, and so gravity helps set this 
scale This is the archetype of a universe dominated by a domain wall, and the domain 

wall is classically forbidden from collapsing to zero radius. However, the universe is of 
hnite volume, ~ .Rmin iiear minimum, so the wall is subject to decay by global quantum 
huctuations, in which the vacuum state in one whole half of the universe jumps to the same 
state as the other half. Clearly this decay process has a hnite, albeit large, barrier factor 
~ . Thus, the domain wall decay problem becomes a problem in quantum gravity. 

The decay process can more particularly be regarded as follows. A second closed do¬ 
main wall nucleates at zero size in the original universe, and the two walls then approach 
each other by quantum tunnelling. When the two walls meet, they annihilate into pure 
held energy. Figure 1 illustrates the 1-wall VIS universe itself, and also the 2-wall decay 
process. For this reason we study the spherically symmetric 2-wall system, coupled to a 
gravitational held, in this paper. An important technical ingredient in this study is the re¬ 
sult of Thiemann and Kastrup [Q, who found an elegant pair of canonical variables (T, M) 
for spherically symmetric gravitational held conhgurations. Here we also hnd compatible 
canonical variables for domain walls. 

Quantum tunneling of domains is already well known in condensed matter physics, and 
has been studied both theoretically and experimentally; see e.g, This gives hope that 

similar processes can be understood in cosmology. 

Section H is devoted to deriving the hrst main technical result of this paper, a simple 
form of the ehective action for the 2-wall system, obtained by integrating out the spherically 
symmetric gravitational held: 
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Nl 

,^R\ - 4\ 

/l — 2M/Ri sin^ — 

- 2M//2i j 

Nl 

I^^Rl - 4\ 

/l — 2M/R 2 sin^ ^jJ2 — — 

- 2M/i?2 j 


Here -Ri, R 2 are the radii of the two walls (dehned by Area= dvri?^), 'ipi and ^’2 are cer¬ 
tain (imaginary) time coordinates along the world sheets of the two walls, (T, M) are the 
Thiemann-Kastrnp variables for the region between the two walls, and /i = dvrcr. Qnantiza- 
tion of this action is straightforward and is carried ont in Sect. Ill; then bonndary conditions 
are set, and the second main resnlt is then derived, the qnantnm tnnnelling probability for 
domain wall decay. 


P ~ exp . 


( 1 ) 


Section IV presents a slightly different qnantization, in 
Section V discnsses the results and compares to paper 
method of Hamiltonian quantization |]9|-|Tl • 


which intrepretation is a bit clearer. 
I. In Appendix A we review Dirac’s 












FIGURES 




FIG. 1. Spaces dominated by one domain wall and by two domain walls. The spaces are closed 
and topologically heavy curves denote closed 5^ domain walls, (a) A space slice of the VIS 
spacetime. Regions Vi and Vq are slices of flat spacetime, meeting at the domain wall, (b) A space 
slice in the quantum tunnelling regime. The second wall has nucleated at zero size and the two 
walls are tunnelling toward toward each other; when they meet they will annihilate into pure field 
energy. Regions Vi, Vq and V 2 are slices of flat Euclidean space. 


II. CANONICAL ACTION FOR THE TWO WALL SYSTEM 


Dirac first applied his method of Hamiltonian quantization (described in Appendix A) to 
general relativity in |]^, and the theory was soon developed in greater detail by Arnowitt, 
Deser and Misner |]^. It is not the purpose of the present study to review the details of 


this subject, however, but rather to use the canonical formalism as an alternate approach 
to the instanton calculation of paper I, which was shown to have a certain pathology. The 
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reader is referred to the literature for a complete treatment of canonical quantum gravity 

0,111 • 

The quantization of a spherically symmetric spacetime with a hnite number of degrees 
of freedom is a minisuperspace model known as the Berger-Chitre-Moncrief-Nutku (BCMN) 
model |T5Hl7| . We will be quantizing the spherically symmetric domain-wall spacetime (VIS 
spacetime p|,||) introduced previously, while allowing for the possibility that a second domain 


wall tunnels from zero size and annihilates with the existing one. Q We therefore expect that 
there will be just two degrees of freedom in the problem, corresponding to the radii of the two 
domain walls. A related problem with one degree of freedom, the canonical quantization of 
a spherical bubble of false-vacuum, was worked out in the WKB approximation by Fischler, 
Morgan and Polchinski [Q. This same problem was earlier studied using the Euclidean 


approach by Blau, Guendelman and Guth |]^, who found a pathology very similar to what 
was encountered in the instanton calculation of our previous paper Q]. 


A. The First Order Action 


The action for gravity plus domain walls is given in the thin-wall limit by 
" = /walls 

where again p/dvr is a, the energy per unit area of a domain wall, and where the sperically 
symmetric metric as 

ds^ = -{NUtf + L^{dr + N^f + R^{de^ + sin^ Odcj)^). (3) 


In what follows, we will take A = 0, corresponding to the case of pure domain walls. Also, 
except where noted, we work in units where G = 1. The canonical coordinates, which are 
functions of (r, t), are (A*, A’’, L, i?, ri, r 2 ), where ri and r 2 are the radial coordinates of the 
two domain walls. Dehning the conjugate momenta as usual, the Hamiltonian form of this 
action is 


where 18 


S = J dt piri+p 2 r 2 + j dr {ttiL +'KrR-N^H t-N''Hr) 


Ht = 


LttI TTlTTr 1 
R ^ 2 


Pi 


2i?2 


2R [R'V R'^ ■ 

t[t) 


1/2 


E 5(G - r) ( ^ + p^R^ 


i=i,2 


L2 


Hr = R'tir - Lvr^ - 6{rj - r)pj. 

J=l,2 


(4) 


(5) 


^One might more generally study non-spherically-symmetric tunnelling configurations. We expect 
these to be no more probable, and to affect the prefactor but not the exponential barrier factor in 
our result. However, we have no proof of these expectations. 
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This action is generally covariant under coordinate transformations of {t,r). 

Since there are no time derivatives of A^t and in the action, the primary constraints 
are 


TTtv* — TIN’’ — 0, 


( 6 ) 


which are hrst-class. However, since the Poisson brackets between these constraints and the 
full Hamiltonian do not vanish, there are the secondary constraints 

nt = nr = 0, (7) 

which are also first-class, and which generate coordinate transformations of {t, r). Assuming 
R{r) to be continuous and 'itl,r to be free of delta-functions at each wall, integration of these 
secondary constraints across a wall implies the following jump conditions at the surface: 

A>r, = -p 

= ( 8 ) 

where E = {p‘^ + evaluated at the wall. 

To implement these constraints in the 1-wall system, Fischler, Morgan and Polchinski 
T8| followed the Dirac approach (see the Appendix A) to hnd a wave function satisfying 

ttat* |T) = TTjv^ Id') = 0, (9) 

Tfi |T) = K 1^) = 0. (10) 

The hrst pair of these relations simply says that the wave function is independent of the 
lapse and the shift; the second pair will generate the dynamics of the wave function. 

As has been mentioned, it is generally true that the hrst-class constraints are in one- 
to-one correspondence with the gauge symmetries of the theory; the existence of four such 
constraints in the present case therefore indicates that there are four gauge degrees of free¬ 
dom. Two of these correspond to the invariance of the theory under different choices of the 
lapse and shift functions; we are also free to hx the time slicing and radial parametrization 
through gauge choices. Whether or not one hxes this part of the gauge before quantization 
distinguishes Dirac quantization from ADM quantization: Dirac’s procedure, involving no 
gauge hxing, leads in principle to the wave function for all possible time slicings and radial 
parametrizations, whereas in the ADM procedure one hxes the gauge before quantization 
and winds up with the wave function only for a given slicing and parametrization. The 
Dirac procedure is generally more unwieldy than the ADM method; however, one must take 
care that possible quantum behavior is not ruled out by a premature gauge choice. Indeed, 
it was shown in |jT^ that overzealous gauge hxing may lead to the inadvertent exclusion of 
parts of the quantum dynamics. 

We will therefore take a hybrid Dirac-ADM approach. Roughly speaking, we will “inte¬ 
grate out the gravitational held”: We will hx the radial coordinate and take a hxed family 
of time slices, and then solve the constraint equations in the three vacuum regions separated 
by the two walls, to reduce the action to an ehective action which exclusively involves wall 
degrees of freedom. Then we will implement Eq. (IT) solely at the walls. 


(In fact, the gauge hxing is a convenience but not a necessity for this problem. A future 
paper in this series will present a [nearly] gauge invariant derivation of the ehective action 
for n walls in spherical symmetry.) 
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B. Gauge Fixing 


We first look for a solution to the constraints as follows. Let the radial coordinate take 
the range 0 < r < rs, and divide the (compact) space into three regions. There are two 
centers of spherical symmetry, located at r = 0 and r = r^, corresponding to the centers of 
the two spherical domain walls, and in addition there is a middle region between the two 
walls, where ri < r < r 2 . The walls themselves are located at r = ri and r = r 2 . We will 
refer to the three regions as region Vi, Vq, V 2 respectively: 


Hi : 

0 < r < ri. 


Id): 

Ti <r < r 2 , 

(11) 

C 2 : 

r 2 < r < To- 



We now £x the radial parametrization everywhere by imposing the coordinate gauge condi¬ 
tion 


L = l, 


( 12 ) 


and then impose the slicing condition 


r Rttji = 27rL, Vl, 

< Rttr = ttl, Vo, 

[ RtIr = 271l, V 2 . 


(13) 


C. Solution of the Constraints 


Combining these conditions with the spatial constraint equation 


leads to the solutions 


Rr 


= 


TTl = 


R’tir - 

ttl 

2ikiR, 

Vi, 

ico, 

K), 

2ik2R, 

1^2, 

ikiR^, 

Id, 

icoR, 

K), 

ik2R‘^, 

Id; 


(14) 


(15) 


where ki, cq, k 2 are constants of integration. Here and below, the constants of integration 
that appear will become our degrees of freedom, and should all be understood as functions 
of time in the dynamical problem. The factors of i have been chosen appropriate to the 
classically forbidden, tunnelling regime. 

The Hamiltonian constraint Rt = ^ then becomes 

r SkjR"^ - 1 + 2RR" + R'"^, Hi, 

{] = 2Rt = \cl-l + 2RR" + Ho, (16) 

\ SkjR^ - 1 + 2RR'' + R'^, H 2 . 
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These equations have the general solutions 


dR 


' - klR^ 

in region Vi (assuming regularity at the origin r = 0 ); 

/■ dR 

r - Vo = 


(17) 


\/l - c§ - 2M/R 

in region Vq, where M and tq are further constants of integration; and 

f dR 


(18) 


r — = 


1 - 


(19) 


in region V 3 (assuming regularity at the anti-origin r = r^) where is another constant of 
integration. 

In region Vi, requiring i?(0) = 0 gives the solution 


Rir) = —sm(kir), 0 < r < ri. 
ki 

Similarly, in region V 2 , requiring i?(r 3 ) = 0 gives the solution 


R{r) = — sin[/i; 2(’"3 ~ ’")]5 r 2 < r < r^. 
K2 


( 20 ) 


( 21 ) 


In region Vq, dehning 


Co = sin 6*0 ( 22 ) 

we leave the solution in implicit form as 

r dR . . 

r - ro = / . ri<r < ra- (23) 

■> ^Jcos'^9o-2M/R 

Here the constants tq and are hxed in terms of the other variables by our requirement 
that R{r) be continuous across each wall. Eqs. (^)- (l^) represent a general solution to 
the Hamiltonian and spatial constraints, parametrized by (0o, ^ 2 , Ta, M). It will be 

convenient in what follows to dehne 

Rj=Rirj) (j = l, 2 ), 

91 = kin, (24) 

02 = vr - ^ 2(^3 - ra); 


here and throughout, the index j = 1, 2 runs over the two walls. We take as our independent 
parameters the set ( 0 o, 6 * 1 , 6 * 2 , Tf, ri, ra). 

Given the above solution to the constraints, it follows from Eq. ([ISD that 
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TTl = 


and 


VTiJ = 


iRsm{kir), 

^ 1 , 


iR sin 00, 

K), 

(25) 

ii?sin[/c2(?’3 — ?")], 

^ 2 , 


2i sin(fcir). 

^ 1 , 


i sin 00 , 

K), 

(26) 

2 isin[fc 2 (r 3 - r)]. 

4"2. 



D. Reduction of the Action 

We have now found a 6-paranieter family of solutions which identically solve the con¬ 
straints Tit^r = 0 everywhere except at the 2 walls, while at the 2 walls we still have 2 
canonical momenta pj, and 4 constraints, the jump conditions (^). Our reduced phase space 
is now of hnite dimension 8, with coordinates (^o, 6*i, 02, ri, r 2 , M,pi,p 2 )- Our system is 
described as a time dependent point in phase space obeying the constraints. 

We can therefore write the total action as 

S = J dt^^pif-i + p2r2 + J dr^niL + tirR - N^Ht - N'^Hr)^ 

= /*{ E (pn - y' [Ei + RiiAR')^] - iVJ [p, + 

S=l,2 

+ J dr{TTLL + TrRR)y (27) 

where Ei and E 2 , obtained from the dehnition Ej = (p'j + are now given by 

f/Rl — Rl{sm 01 — sin 0o 


E, = 

E2 = 


/i^i ?2 ~ 7?2(sin 00 — sin 02 )^ 


1/2 


11/2 


(28) 


and where, using Eqs. 


(A7?')i = cos 01 — Y^cos^ 00 — 2M/Ri, 

{A71r)i = —-ii?i(sin0i — sin0o), 

{AR ')2 = — cos 02 -|- cos^ 00 — 2M/ R 2 , 
{Attl )2 = -ii? 2 (sin 0 o - sin 02 ). 

Next we calculate the gravitational contribution to the action, i.e., 


(29) 


(30) 


In the gauge L = 1, we can neglect the hrst term. The second term is calculated as follows: 
write Eqs. (|^)-(^) in the form 
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R(r) = AeA + BeB + cec, 


( 31 ) 


where A, B, C are just different names for R{r) in the three regions, and 

©A = 0(n - r), 

©B = 0(r - ri) - 0(r - r 2 ), (32) 

©c = 0(r -rs) 

are step functions for the three regions. Similarly, write Eq. (PB[ ) as 

ttr = 2ikiAQA + i sin 6*o0 b + 2ik2CQc- (33) 

Then it follows that 


ttrR = i/ci(ff^)*0A + i sin OqBQb + ik 2 {C‘^)*Qc- (34) 

Integrating by parts and using the fact that 

©A = ri5{ri - r) 

©s =-ri5(r - ri) + r2<5(r - r2) (35) 

©c = - r2). 


one finds 

Sg = J dt^^ri{—ikiRl + isinOoRi) + r2{ik2Rl — ism9oR2) 

+ J (ir(——-icos^o^o-B©^ — *^ 2 C'^©c)| (36) 

= (A7rL)iri + (A7rL)2r2 + J dr{-ikiA‘^QA - icosOoOoBQB - ik2C‘^&c)^- 

where 9i = fciri, 6*2 = vr —A) 2(’"3 —'^ 2 )- Combining Eq. (^) with Eq. (^Tf) , the action becomes 


where 


S' = y dtdr{—ikiA?QA — icos9o9oBQB — ik2C^Qc) 



(37) 


Pi=Pi + (A 7 rL)i= Pi - iRi{sm9i - sin^o), 

P 2 =P 2 + (A 7 rL) 2 = P 2 - iR2{An9o - sin 6 ^ 2 ). (38) 


Now consider the terms involving A and C in Eq. ( 0 - Performing the radial integration 
leads to 


J dtdri-iklA^QA - ikC^Qc) 
= [ dt 


z/l/2 

7 ^(ri - Ri cos^i) - 7 ^(r ’3 - r 2 + i ?2 cos 6 * 2 ) 


(39) 
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Furthermore, using the dehnitions kiRi = sin6*i, 6 i = kiVi and /C2-R2 = sin6*2, 6*2 = 
TT — /c2(^’3 ~ 1^2), one can show that 


(ri - Ri cos ^i) = Rl 6 i + 

dt 

— r 2 + R 2 cos 62 ) = —R ^92 + 

where 

R2 

“ = - 0 - 2/1 (^1 - cos^i), 

z sin ui 

r2 

■ In - 6^2 + sin 02 COS 02 ) . 

2 sin 02 

Hence we can write, dropping the total time derivatives, 

J dtdr{—iklA‘^QA — ik2C^Qc) = J dt{—iR\0i/2 + iR^92/2). 


The term involving B can be rewritten as 


dtdr{—i cos 9 q9oBQb) = / dt{—i cos 9o9o) 


r2 RdR 


R' 


= / dt{-i9o) 




(40) 


(41) 


(42) 


(43) 


where 


F{M,Ri,R 2 ) 



cos 00 

^/co^9^^^2M/R 


(44) 


Combining Eqs. (|37[), (^2]) and (|43|), we hnd the reduced action 

S = J dti^iRl9i/2 + iRl92/2 - 19^ (rII2 - R\l2 + F(M, i?i, R^^ 
- E {N^AEj+R,my + 


f-i - n: 


Pj 


i=i,2 


}}■ 


(45) 


Consider the last terms inside the sum over walls, the terms proportional to pj, in Eq. (^51). 
From the spatial constraints, pj = 0 for all time, so these terms make no contribution to 
the equations of motion, and can be dropped from the action. The old, gauge-dependent 
canonical wall coordinates rj therefore disappear from the action, in favor of the gauge 
invariant wall quantities Rj = R{rj). Moreover, the old, gauge-dependent wall canonical 
momenta pj have completely decoupled from the remainder of the action, of their own accord, 
and can be dropped henceforth, along with the spatial constraints. The quantities i9i, i 02 
now act as canonical coordinates for the two walls, and the quantities R‘j/2 = R?{rj)/2 now 
act as gauge-invariant canonically conjugate momenta for the two walls. The quantity i9o is 
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a single remaining canonical coordinate of the gravitational field in region Vq, and M is some 
kind of momentum belonging to it, though not canonically conjugate. Evidently, spherically 
symmetric gravity is one of those simple gauge theories wherein the unphysical degrees of 
freedom decouple of their own accord when appropriate canonical coordinates are chosen. 

The Hamiltonian constraints at the walls can also be simplihed. They now read 




f/Ri — Rl{sm9i — 
~ Rl{sm9o — 


sin 9oY 
sin 92 ^ 


+ i?i(Y^cos^'0o'--'2iw7^ ~ cos^i)! 

+ i? 2 (cos 6'2 - \J cos^ 9o - 2M/R 2 )|. 


(46) 


(47) 


We have here the usual awkwardness that the Ej from Eqs. (^81) contain a square root, and 
in order to obtain a simpler quantum mechanical system, we take the usual remedy and 
“square out” these constraints as follows. Dehning 

TV] ^ RfN^^ [E, - R,iAR%] {j = 1, 2) (48) 


the constraints become 





eVr^ar!) 


-Nl 

-Nl 


fX^Rl — 1 + 2(cos 9i \J cos^ 6*0 — 2M/Rq + sin 6*1 sin 6 * 0 ) 
~ 1 + 2(cos 6*2 \J cos^ 6*0 — 2M/i?o + sin 6*2 sin 6 * 0 ) 


1 + 2 M/i?i 

1 + 2M/i?2 


(49) 


The extra factors introduced into the constraint by Eq. (|^) never vanish in the classical 
regime, and make no difference to the classical equations of motion. Moreover they never 
vanish in the quantum tunnelling regime and so at most affect the prefactor in the tunnelling 
calculation. The effective action is now 




jdt 

-Nl 

-Nl 


iRl9i/2 + iRl92/2 - (rI/ 2 - -Rl/2 + F(M, i?i, R 2 )) 

fi^Rl — 1 + 2(cos 9i cos^ 9q — 2M/Rq + sin 9i sin 6 ^ 0 ) — 1 + 2M/Ri 

/r^i ?2 ~ 1 + 2 (cos 02 cos^ 9o — 2 M/Rq + sin 6^2 sin 6 ^ 0 ) — 1 + 2 M/R 2 


(50) 


and the reduced phase space is now 6 -dimensional with coordinates ( 6 * 0 , 9 i,92, M, Ri, R 2 ). 


E. A Cyclic Time Coordinate 

We now wish to hnd a further coordinate transformation in phase space to canonical 
form. In Eq. (|5(]D, the quantity M which appears in the integral expression for E{Ri, R 2 ) is 
the Schwarzschild mass of the region of spacetime between the two walls. Hence we expect, 
and conhrm, that variations of this action lead to M = 0 as the equation of motion for M. 
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However, it is clear that there is no cyclic coordinate present for which M is the conjugate 
momentum. M is part of the conjugate momentum of but 6q is manifestly not cyclic in 
the action. It is thus desirable to hnd the coordinate to which M is canonically conjugate, 
and which would therefore be cyclic in the action. 

Recently, Thiemann and Kastrup found ^ that such a canonically conjugate pair of 
observables can always be found for spherically symmetric held conhgurations; see also 
PIU^^. They worked in the Ashtekar approach to canonical quantum gravity, but their 
result is general. They showed that for the line element 


ds^ = -{N^dtf + L^{dr + + R^de^ + sin^ Odcj)'^), 

the time variable conjugate to the Schwarzschild mass is given by 

T = - J dr{l - 2M/R)-^^{dR/drf - L‘^(l - 2M/R), 

which in our setting is given by 

dR 


(51) 


(62) 


T = —i sin 9n 


rR2 


-^cos^ 00 - 2M/R(1 - 2M/R) 

We would like to relate this to the function F dehned by Eq. (|4^ . To do so, note that 

1 


(53) 


dT . 

— = -i sm 00 

dRi,2 


(l-f)(cos2 0o-f)V2 
dR 


R2 


J Ri 


dT f 

rin * COS 00 / , 

900 J (cOs2 00 - ^)3/2 


(54) 


dT 


= —zsin0( 


R 

dR/R 


+2 


° ' (l-^)(cos 2 0o-^)3/2 

dR/R 


(1_^)2(cos2 0 o-^) 1/2 


We can integrate by parts in the hrst term in dT /9M, after writing the integrand as 
[R/{1 - 2M/R)] ■ [dR/[R2(cos2 0o - 2M/Rf/^] to give 


dT 


i sin 9oRj 


M(l-f^)(cos2 0o-f^)V2 


R2 


Ri 


T 

^ M' 


( 55 ) 


The time derivative of T is then given by 


dT ■ dT ■ dT ■ dT ■ 


( 66 ) 


Furthermore, we can integrate F by parts to hnd 
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RdR 




+ 




-|H 2 

I 

-IKi 


i9o cos 9qM 


dR 

(cos2 0o-^)3/2 


i9o cos 9oRj 


2 (cos^ 6*0 


2M^l/2 


Ri 


R 2 


+ 


'^1^0 


J Ri 


we find that T and F are related by 


Combining Eqs. 

^{MT-MT) = -i9oF{Ri,R2,cos9o,M) - 


2 

+ 




1^2 

I 

JRi 

R2 

Ri 


2^Jco^~9o^^2Mj^j 


^R] cos 00^0 + sin 9 ^^^— 

\ ^ ^ ^ “ 1 - 2M/Rj 


R 2 


-* Ri 


( 57 ) 


(58) 


Up to a total time derivative, therefore, we jnst get the desired term MT in the action, pins 
some messy terms at the two walls. 


F. Canonical Wall Coordinates 


Our next task is to hnd compatible canonical coordinates {'ipj, Rj) for each wall. Consider 
the N\ constraint at wall 1. The angles 9o and 6 *i both enter into this constraint, but we 
would like to reduce this to a dependence upon a single angle 'ipi. Guessing the answer, we 
can do this by dehning 'ipi by 

sin( 6 *i + 2%jji) = sin 6 *o/\/l — 2M/Ri (59) 

which implies 

cos( 6 *i + 2^jJl) = sj cos^ 6*0 — 2M/Ri /— 2M/Ri (60) 

Then the constraint becomes 

Rj = - 2M/Ri sin^ + (1 - ^Jl - 2M/Fl)^ (61) 

9o and 6 *i themselves having disappeared as desired. Also, by time-differentiating the deh- 
nition of ipi we hnd 


9\ -\- 2'ij)i 


1 

Y^cos^ 9o - 2M/Ri 


^cos 6 * 0^0 + sin 6 *o 


MRi-MRi \ 
Rl{l-2M/R,)) 


(62) 


an expression which furnishes exactly the wall terms that are needed in Eq. (|58|)for T. 
Similar equations hold at wall 2: 


sin( 6*2 - 2 ip 2 ) = sin 6 *o/- 2 M/R 2 , 


(63) 
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which implies 


(64) 


cos(92 - 2 i/> 2 ) = ■Jcos‘e„ - 2M/R2/-J1 - 2M/R2. 
The wall 2 constraint is 

= 4-^ 1 — 2M / R 2 sin^ '^2 + (1 ~ — 2Tf/i?2)^; 

and for for the time derivative 


62 — 2'^2 


1 

^COs2 00 - 2 M/R 2 


cos + sin 9q 


MR2-MR2 \ 
RI{1-2M/R2)) 


(65) 


( 66 ) 


G. The Effective Action for Two Domain Walls 

Pntting all this into onr previons action gives 

S = J dtf^iRliji + iRlij2 + Mf 


Nl 

- 4\ 

/l — 2M/Ri sin^ -01 — — 

- 2M/i?i j 

. 9l 

Nl 

I^^Rl - 4\ 

/l — 2 M/R 2 sin^ '02 — (^1 — 

— 2M/ i?2 ] 


as the effective action for two domain walls, with canonically conjngate phase space coor¬ 
dinates 0 (T,i'ipi,i'ilj 2 , M, Rf, R 2 ), in which the gravitational held is rednced to the single 
degree of freedom T, and in which there is one degree of freedom 'ijjj belonging to each of 
the two walls. There remain two constraints. 

We note the following simplihcations that have taken place along the way in the derivation 
of this action; 1) The original, non-gange-invariant, canonical coordinates {rj,Pj) of the walls 
have dropped ont completely, in favor of effective, gange-invariant, canonical coordinates 
[Rj^Rj). 2) The gravitational degrees of freedom have all been integrated ont, to leave 
behind jnst (T, M). 3) The gange dependent coordinates iOj which arose dnring solntion of 
the constraints have likewise disappeared. 


^ In this action, the coordinates ipj and factors of i have been chosen appropriate to a tun¬ 
nelling problem. The turning points to the classically allowed regimes of phase space are at 
'ipj = ±7r/2, ±37r/2, ±57r/2... and, if desired, the coordinates can be analytically continued at 
those points. See Sect. IV. 
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III. QUANTIZATION AND THE WAVE FUNCTION OF THE TWO WALL 

SYSTEM 

Starting from the action given by Eq. (|^, we now proceed to quantize our 2-wall 
system, set boundary conditions appropriate to the initial state of the system, solve for 
the 2-wall wave function, '^ 2 , E) and so determine the decay amplitude. This system 

can be quantized exactly and it is unnecessary to resort to the WKB approximation.^ The 
action, Eq. (0. contains two constraints, which can be treated by Dirac quantization (See 
Appendix A). 


A. Reduction to M = 0 


In the initial state — a single VIS domain wall at minimum radius — we have M=0. 
Thus the boundary value of the wave function T is independent of T, since (T, M) are 
canonically conjugate P,pO|- p^] . But M commutes with the constraints since T is a cyclic 
coordinate, so T must be independent of T throughout: 

= T('0i,V’2) (68) 


Thus the domain wall decay problem reduces to a 4-dimensional phase space in 
(- 01 , ' 02 , -R 15 -R 2 ) with action 


S = I dt^ iRf'ipi + zi?202 — Nl — dsin"* 0i 


m 


“ 4sin^ 02 


The canonical coordinates are now (z0i,z02) and the canonical momenta are 

TTi = Rl, 

712 = Rl, 


SO that the constraints can be rewritten 


H] = —if/Tii — 4sin^(0i) = 0, 

— 4sin^(02) = 0. 


(69) 


(70) 


(71) 


B. Quantization 

Quantizing by Dirac’s procedure, we promote the constraints to operator equations to 
define the physical state space. 


^But we have indeed made some choices in the formulation of the constraints, which make no 
difference to the classical system, but which do affect details the quantum system beyond the 
WKB approximation. See Sect. IV for further elaboration. 
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H] 14') = 0, 

W? 14-) = 0. 

We thus take momenta as differential operators 


( 72 ) 


VTl 


. d 


d 


diipi d-ipi ’ 

. d d 


^ di'ilj2 ' dilj2 ’ 

acting on the wave function -^ 2 ), and we arrive at the two wave equations 

+ 4sin2(^7i) d'(V^i,V’2) = 0, 

oi>i J 

+ 4sin2(V>2) = 0, 

di>2 


(73) 


(74) 


Since these equations are uncoupled, the wave function separates as the product of 1-wall 
wave functions, i.e., 


4'('/’i,V’2) ^ 4'i(i/ii)4'2(V’2)- 


(75) 


It may seem surprising that the two walls are uncoupled from each other; however this is a 
necessary consequence of Birkhoff’s theorem and conservation of M. We point out two facts: 
First, the two walls are topologically coupled by residing in the same space; second, there 
will be a local effective coupling between the two walls (he., a ^-function coupling) which 
will mediate their annihilation into pure field energy, that we have missed in the thin-wall 
approximation. Then our wave equations become 




4'(</>,) = 0. 


(76) 


C. Boundary Conditions 

We now set appropriate boundary conditions for the wave functions whose dynamics 
are given exactly by Eq. dlD- These equations take the slightly unusual form of first- 
order equations, not the usual second-order Schrodinger equation. First-order equations 
can, however, be routinely handled by Dirac quantization |p. Q It is clear that, in the 
tunnelling regime, solutions will die exponentially as ipj increases. The usual freedom to 


slightly different, second-order form of the action will be presented in Sect. IV. There is also 
still another way to proceed that also leads to second-order Schrodinger equations. We can put ± 
signs on the radicals in Eq.(^), to impose two separate constraints at each wall instead of one. 
Classically, a solution would be such that one of these constraints vanishes; the wave function, 
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choose either exponentially dying or exponentially growing solutions appears in a slightly 
different way here: We can choose a 'ipj in the initial state corresponding to the desired value 
of Rj, and then evolve either to increasing 'ijjj or to decreasing respectively. Since our 
problem involves tunnelling, we wish to evolve both -0^ to increasing values. 

We now restore the factors of G, where 




/^gut 

V >1 


6 


(78) 


The particular problem we wish to solve involves a pre-existing wall at the turning point 
Ri = 2//iG, and a second wall just nucleating at zero size, = 0. Classically from Eq. (|6^, 

^i = ^sinV^j, (79) 

so translating the Rj to canonical wall coordinates. 


= 7r/2, 

'11^2 = 0, (80) 


which give the initial conditions on the wave function. 


'I'(7r/2,0) ^ Cnuc 


(81) 


where Cnuc is the amplitude to nucleate a zero-size wall, which we take to be ~ 1. 

The tunnelling process then proceeds with both 'ijjj increasing, until the walls meet and 
annihilate. This necessitates for the hnal state 


Ri — R 2 , 

or (taking into account the multiple-valuedness of sin) 


V’l = 


^2 

vr - -02 


The allowed domain for the canonical coordinates is 

7r/2 < "01 < TT, 

0 < ipi < 1^12 


(82) 


(83) 


(84) 


representing a linear combination of such allowed states, will only be annihilated by the product 
of the constraints. Hence for the purpose of quantization, we write 

n{ = (77) 

This eventually leads to equivalent results (up to factor ordering ambiguities) after careful impo¬ 
sition of boundary conditions. 
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so the final condition must be 


'lpl=TT- '02- (85) 

The final state is not a turning point of the classical action, Eq. (|6^), because the annihilation 
process is missed in this action, as mentioned above. Presumably a small-scale turning point 
would appear if we worked beyond the thin-wall approximation, but we will not pursue this 
point. Figure 2 displays the configuration space of the 2-wall system in the tunnelling regime. 



FIG. 2. The configuration space of the two wall system in the tunnelling regime. The coordi¬ 
nates are ipi and 'ijj 2 ■ Initial conditions are shown as • at the lower left. The diagonal dotted line 
represents the locus where the walls collide and annihiliate into pure field energy. 


D. Solution of the Wave Equations 

The wave equations (^,^) comprise two independent first-order equations in the (V’l, '?/’ 2 ) 
plane, which are well posed and have a unique solution under the the initial conditions (|^ 
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on the domain (|^). This solution, for the wave function of the 2-wall system, is then 


= Cnucexp 




2 ijji + 2 ip 2 — vr — sin 2'^i — sin 2 ip 2 


( 86 ) 


Evaluating this wave function under the hnal condition, Eq. (|85[) gives 

'hfinai('0i, vr - -01) = Cnucexp (-7r//r^G^) (87) 

for the hnal state wave function. The corresponding probability of tunnelling to the hnal 
state is 


T^finai = |cnuc|^ exp (-27r//i^G^) (88) 

The most remarkable feature of the result of this paper, Eq. (^81) , is its independence of the 
hnal value of "01 or 'ip 2 - This means that the two walls may collide and annihilate at any 
value of the hnal radius in the kinematically allowed range, 

0 < i?final < 2//iG, (89) 

with equal probability. At hrst this may seem surprising, but we argue that it is as expected. 

Consider the following toy problem. A particle and an antiparticle move in a potential 
that is identical for both particles. (For instance, a proton and an antiproton move in a 
gravitational potential.) There is a potential barrier present, and in the initial state, the two 
particles are on opposite sides of this barrier. They may tunnel toward each other through 
the barrier, and annihilate if they meet. The question now is, what is the most probable 
location for the annihilation? We encourge the reader to stop reading at this point, guess 
the answer, and then work it out. 

The answer is that annihilation is equally probable at any location within the barrier, and 
the annihilation probability is just given by the total barrier factor for single-particle pene¬ 
tration. We argue that annihilation of the 2-wall system is no different, justifying our result. 
However, this system is not easy to interpret in the canonical variables Therefore 

we will also give a slightly different quantization for the 1-all system, with application to 
the 2-wall system, after dehning some new canonical variables. 


IV. A FURTHER METHOD OF QUANTIZATION 

A better pair of variables (Q, P) for the 1-walll system can be obtained by dehning 

X = i{T^/‘2-'4)) (90) 

and then carring out the following canonical transformation: 

Q = \/2R cosh X, 

P = v^i^sinhy, (91) 

which entail 
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iR^ip = PQ + (total time derivative), 


R — AR sin xjj = ^ R — AR cosh X; 

= ^{Q^-P^f-2Q^, 

V2, 


V2, 


(We are again setting G = 1.) The 1-wall action then can be rewritten as 

S= fdt\iR^-N^ fi^R^ - aJi - 2M/Ri sin^ (i-Ji- 2M/Ri 


I dt IpQ - N‘ 11(152 - p2) _ 21q\\= ldt{PQ- N‘H,} , 


( 92 ) 


(93) 


where the constraint is redehned by 


nt = ntR-^ + ’ 

iV* = iV*/?-2 |^i(Q2_p2^ + ^Q^ . (94) 

Once more, the extra factors introduced into the constraint by Eq. (|^) never vanish in the 
classical regime, and make no difference to the classical equations of motion; and they never 
vanish in the quantum tunnelling regime and so at most affect the prefactor in the tunnelling 
calculation. The phase space (Q, P) of the 1-wall action in the form (|^) can be described 
as follows. The classically allowed regime is 

Q real, 0 < Q < cx); P real, 0 < |P| < Q; (95) 

while the classically forbidden, or quantum tunnelling, regime is 

Q real, — oo < Q < oo] P imaginary, — ioo < P < ioo. (96) 


The two regimes meet at P = 0. 

The action (p3| ) is now entirely straightforward to quantize as a 1-dimensional particle 
system. We take the 1-wall wave function as T(Q) and use 


P 



( 97 ) 


to write the constraint as a Schrodinger equation TijT = 0 or 

where the potential is an “upside-down harmonic oscillator”. 
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( 99 ) 


V{Q) 


-Q^ + —<3, 



Note that the constraint says that ^ mnst be the zero eigenfnnction of Tit] the rest of 
the spectrnm of Tit is not present. Eqnation (|^) is hypergeo metric and its solntions are 
parabolic cylinder fnnctions, bnt we will not pursne the details. Fignre 3 shows the potential 
V(Q) and the resnltant dynamics. 



FIG. 3. Dynamics of the quantized VIS solution. The domain wall can be viewed as a particle 
moving in one dimension Q, under the influence of a potential V{Q) (curve). The energy is 
constrained to be 0. The turning point is at Q = 2\/2//r; to its right is the classically allowed 
regime, and to its left is the classically forbidden, or quantum tunnelling, regime. The dashed 
horizontal line represents the classical motion of the VIS wall, and the dotted line represents 
quantum tunnelling. 
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Tunnelling from Q = 0 to Q = 2a/2/ /i is the “creation of the VIS universe from nothing”, 
while tunnelling in the opposite direction is “annihilation of the VIS universe into nothing” 


(c/. @). 


To study the quantum decay of the VIS universe into a universe containing pure field 
energy, our subject in this paper, two copies of the 1-waIl system must be coupled to make 
the 2-wall system with configuration space {Qi,Q 2 )- Boundary conditions are as follows: 
For wall 1, impose along the Qi-axis purely right-going boundary conditions on the left. 
For wall 2, impose along the Q 2 -Sixis purely left-going boundary conditions on the right. 
Annihilation can occur at any Qi = Q 2 between the endpoints 0 and 2\/2/jjL. The WKB 
approximation gives the same exponential barrier factor P ~ exp(—7r//r^G^) as in Eq. (|88D 
above, and the same main result appears, that the decay probability is independent of wall 
radius at annihilation (or final Q). We will not pursue the details beyond the WKB regime; 
presumably the prefactor will differ. 


V. CONCLUSION 


The main results of this paper are: 

1. The 1-wall VIS universe does undergo quantum decay into a universe containing pure 
held energy, with some small probability. 

2. The decay process can be treated as the nucleation of a second domain wall at zero size, 
followed by quantum tunnelling of the two walls toward each other, and anniliation 
when they meet. 

3. The 2-wall system can be treated in a Hamiltonian approach, using a simple action in 
the 2-wall phase space, and Dirac quantization. 

4. The decay probability for the VIS universe is independent of the radius of the hnal 
universe (up to prefactors), and is given by Equation ([881) . 

Conclusions 1 and 2 agree with paper I |l[], which employs an instanton approach to 
replace 3. However, conclusion 4 differs strongly from paper I, which predicted a unique 
value 




( 100 ) 


for the annihilation radius. Furthermore, the probability differs: The above value, Eq. (|88D 
differs from the value of paper I, Eq. (1.88), which is 



( 101 ) 


We interpret this disagreement as an incorrect result of the instanton approach to this 
problem (at least as done in paper I). 

In fact, the two results can be reconciled, if we imitate the Hamiltonian calculation and 
flip some signs in the instanton calculation, in a way that seems ad hoc in the instanton con¬ 
text. In particular, if the four segments in the n = 2 instanton are weighted (-1-1, —1, -|-1, —1) 
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in calculating the action (see paper I), rather than using the (+1, +1, +1, +1) that was pre¬ 
viously motivated both by the standard methodology of Euclidean quantum gravity, and by 
the 2-sheeted manifold rule of Farhi, Guth and Guven [^. These sign flips are now moti¬ 
vated by a careful consideration of exponentially growing versus exponentially dying wave 
functions. After the signs are flipped, the instanton decay probability agrees with Eq.(^). 
We conclude that the instanton method as utilized in paper I makes incorrect choices for 
these wave functions. We would like to be able to propose a “modihed rule” for sign weights 
in the instanton calculation that would repair this defect, but have been unable to hnd a 
convincing formuation. 

We leave the correct instanton treatment of the quantum decay of domain walls as an 
open problem. 

Our result ( |88[) also disagrees with the answer one obtains by assuming that the tunnelling 
probability is of the form 


P 


ex 


vih - h), 


( 102 ) 


where Ij and Jj are the Euclidean actions of the instantons which mediate the creation from 
nothing of the hnal state and the initial state, respectively, as is often done. But, our answer 
does happen to be the same as the probability exp(—Jj ) for creation of the initial state alone 
from nothing, as calculated by using the Euclidean VIS solution as the instanton |^. It is 
not clear why the various methods do not agree. 

The barrier factor we have calculated is a function of the dimensionless parameter 
~ (’^gut/'^pl)^- ^ typical theory, ~ 10^“^—10^® GeV; hence is 

expected to be extremely small in most phenomenologically viable models of microphysics. 
However, improbable events can be important in early-universe cosmology, if they lead to a 
universe resembling our own. 

The new universe created by the decay does not yet resemble our own, however. To do 
so it must hrst expand greatly, and then it must homogenize itself. Whether it does so will 
be the subject of a future paper in this series. 
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APPENDIX A: QUANTIZATION OF GAUGE THEORIES 


Dirac hrst worked out the theory of quantizing constrained systems in general and 
general relativity in particular [^, and his pioneering work continues to serve as the foun¬ 
dation of current efforts to canonically quantize gravity. What follows will be a very brief 
review of the main elements of such a quantization scheme, sufficient for the purposes of the 
current study. Many more extensive studies of the subject can be found in the literature; 
see, for example, |[TT|J^J^ . 
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Consider a mechanical system with n degrees of freedom, which is described by the 
Lagrangian L{qi,qi), and where the qi(t), i = are generalized coordinates. The 

canonical momenta are dehned by 


Pi = 


dL 

oqi 


(Al) 


To pnt this action into Hamiltonian form, one seeks to eliminate the velocities in favor 
of the momenta Pi throngh the nse of Eqs. (|Al|). However, in the event that the Hessian 
matrix of L with respect to the velocities has zero determinant, i.e., 


DetHij 


Det 


dqidqj 


0 , 


(A2) 


then not all of the qi can be eliminated in this manner. (This will occnr, for example, if the 
action is linear in one or more of the velocities, and qnadratic in the rest.) In fact, one can 
eliminate exactly R of the g*, where R < n is the rank of l-Lij. After doing so, one is left 
with a set of {n — R) constraint eqnations of the form 


Ca=Pa- fa{qi,Pi) = 0, a = 1,... ,n - R, 


(A3) 


which are known as the primary constraints of the theory. 
The canonical Hamiltonian, 


H, = f^p,qi-L (A4) 

i=l 

is not nniqne on the fnll phase space {qi,Pi), and so one dehnes the new Hamiltonian 

H = H^ + (A5) 

where the A“ are arbitrary fnnctions, or Lagrange mnltipliers. 

Introdnce the notation {u, n} as the Poisson bracket of the fnnctions u(p, q) and n(p, g), 
and let {u^v}' denote a Poisson bracket to which the constraints have been applied after 
the calcnlation of the bracket. Then one divides the primary constraints into two classes, 
according to the algebra of their Poisson brackets. Those constraints whose Poisson bracket 
algebra closes, i.e., for which 


{C^,Cy} = flpC„ (A6) 

or 

{C„,C^r = 0, (A7) 

are known as first-class constraints, and all other constraints are known as second-class. We 
will denote a second-class constraint with a Latin index, e.g., Ca. An important fact is that 
the for each hrst-class constraint there is a corresponding gange symmetry of the theory. 
Since the constraints shonld hold at all times, we reqnire that 
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C'„ = {C'„,if}' = 0, 


(A8) 


which may in some cases lead to inconsistent equations, or which may lead to new relations 
between the phase space variables. In the latter case, these new relations are known as sec¬ 
ondary constraints, which are then classihed as hrst-class or second-class as described above. 
The consistency conditions ( |A8|) must then be checked again, and the process repeated until 
all constraints have been found. 

In order to quantize the theory, we would like to replace the Poisson bracket in the 
classical relations by —i/h times the commutator of the corresponding quantum operators, 
and then impose the constraints as conditions on the state vectors. However, note that 

Call/') = 0. Cli\'<l}) = 0 

^IC„C^]M = 0, (A9) 

which corresponds to the classical relation 

{Ca,Cyy = o. (Aio) 

Hence all of the constraints should be hrst-class in order for the quantization to go through 
in a straightforward way. 

The prescription for eliminating the second-class constraints is as follows. The Dirac 
bracket is dehned by 


{A, B}* = {A, B} - {A, Ca}T-^{C,, B} (All) 

where the algebra of the second-class constraints is 

= (ai2) 

This is a projection of the Poisson bracket onto the second-class constraint surface; therefore 
if the Poisson brackets are replaced with Dirac brackets in the classical analysis, we can 
consistently take Ca = 0 to hold as operator equations in the quantum theory. Finally the 
quantization can proceed, now with the Dirac bracket taking the role of the Poisson bracket 
in the classical theory: 


h d 
Pi 

i dqi 

Ca{Pi,qi)\y) = 0 . 


(A13) 
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